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Abstract 

A new formulation for the Ffowcs Williams- 
Hawkings quadrupole source, which is valid for a 
far-field in-plane observer, is presented. The far-field 
approximation is new and unique in that no further 
approximation of the quadrupole source strength is 
made and integrands with r“^ and r~^ dependence 
are retained. This paper focuses on the develop- 
ment of a retarded-time formulation in which time 
derivatives are analytically taken inside the integrals 
to avoid unnecessary computational work when the 
observer moves with the rotor. The new quadrupole 
formulation is similar to Farassat’s thickness and 
loading formulation lA. Quadrupole noise predic- 
tion is carried out in two parts: a preprocessing stage 
in which the previously computed flow field is inte- 
grated in the direction normal to the rotor disk, and 
a noise computation stage in which quadrupole sur- 
face integrals are evaluated for a particular observer 
position. Preliminary predictions for hover and for- 
ward flight agree well with experimental data. The 
method is robust and requires computer resources 
comparable to thickness and loading noise predic- 
tion. 

Notation 

□ “ wave operator 

c sound speed in undisturbed medium 

dQ element of the collapsing-sphere surface 

dS element of the rotor-blade surface 

f = 0 function that describes the rotor-blade 
surface 

/+ =; 0 surface described by union of rotor-blade 
surface and rotor disk 
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p = 0 surface that describes the collapsing sphere, 
g = T — t + rfc 
H{f) Heaviside function 
li local force that acts on body 

M local velocity vector of source normalized by 
c, with components M,- 
Mi dMi/dr 

Mr Mach number of source in radiation 
direction, M,ft 
Mr MiVi 

Mr Mih 

Mat advancing tip Mach number 

Mh hover tip Mach number 

n unit outward normal vector to surface, with 

components hi 

p' acoustic pressure, p — Po 

Qij quadrupole surface source tensor, symmetric 

Qm M Q ij Aft Afj 

Qmt QijMiVj 

Qmv Qij^i^j 

Qrr Qij^i^j 

Q Mr Qij Af j f’j 

Qrr Qij^i^j 

Qrr Qij^i^j 

r distance between observer and source, 

r = |x - yl 

r unit vector in the radiation direction, with 

components Vi 
R rotor radii 

t observer time 

Tij Lighthill stress tensor, puiUj + {p' — 

(inviscid form) 

Vn local normal velocity of source surface 
X observer position vector, with components x, 

y source position vector, with components y, 

Greek symbols; 

F intersection of collapsing sphere g = 0, and 

source surface f = 0 
S{f) Dirac delta function 

0 angle between h and f 

po density of undisturbed medium 
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r source time 

Subscript: 

ret quantity is evaluated at the retarded time, 

T = t — r/c 

Introduction 

High-speed impulsive (HSI) noise is a particularly 
intense and annoying noise generated by helicopter 
rotors in high-speed forward flight. This HSI noise is 
closely associated with the appearance of shocks and 
transonic flow around the advancing rotor blades. 
The quadrupole sources in the Ffowcs Williams- 
Hawkings (FW-H) equation^ account for nonlinear- 
ities in the vicinity of the rotor blade. These non- 
linearities are of two types, which are described by 

o 

Lighthill."^’ First, the local speed of sound is not 
constant but varies due to particle acceleration. Sec- 
ond, the finite particle velocity near the blade influ- 
ences the velocity of sound propagation. By inclu- 
sion of the quadrupole source, the correct physics 
is mathematically simulated in the acoustic anal- 
ogy. The quadrupole source in the FW-H equation 
was identified by Yu et al."^ as a significant contrib- 
utor to helicopter HSI noise. Hanson and Fink also 
included the quadrupole source for high-speed pro- 
peller noise prediction but found that it was not a 
significant noise source in that application. Even 
though this early work demonstrated the importance 
of the FW-H quadrupole, it has not been routinely 
included in rotor noise predictions because of the 
difficulty in predicting the source strength of the 
Lighthill stress tensor Tij and the lack of a com- 
putationally efficient algorithm for computing the 
quadrupole noise. 

In the past few years, the computation of the tran- 
sonic aerodynamic field around rotor blades has be- 
come feasible; hence, renewed interest in prediction 
of HSI noise has emerged. Yu et al.'^ were the first 
to successfully utilize advances in CFD by approx- 
imating the quadrupole source strength and inte- 
grating in the direction normal to the rotor plane. 
The integration in the normal direction of the ap- 
proximate quadrupole source, which is valid in the 
far field ahead of the helicopter, effectively trans- 
forms the volume integration of the quadrupole into 
a surface integration. More recently, Schultz and 
Splettstoesser,® Schultz et al.,^ and lanniello and 

Q 

De Bernardis have also used this technique with 
good results. Prieur^ and Prieur et al.^® have devel- 
oped a frequency domain method for computing the 
quadrupole noise of hovering rotors that has yielded 
good results. 

Some attempts have been made to numerically 
integrate the entire volume around the blade, ® 


but the computations generally require computer 
resources comparable to those required by un- 
steady three-dimensional computational fluid dy- 
namics (CFD) calculations — significantly more than 
that required for thickness and loading noise pre- 
dictions. Farassat^^ and his colleagues^^’ also 
tried to reduce the computational effort required 
in computing HSI noise; they recognized that the 
appearance of a shock wave coincides with the on- 
set of HSI noise. By assuming that the shock is 
the dominant contributor of quadrupole noise, the 
acoustic sources are mathematically confined to the 
shock surface. When the shock-noise theory was 
implemented, the conclusion that the shock noise 
was a dominant component of the quadrupole source 
was verified.^’’ Nevertheless, the difficulty in accu- 
rately extracting the shock geometry, location, and 
strength from CFD solutions, together with the fact 
that the shock noise alone did not sufficiently charac- 
terize the total quadrupole source contribution, has 
postponed the complete implementation of the the- 
ory. 

The goal of this work is to utilize the far-field 
approximation to the FW-H quadrupole given by 
Brentner and Holland^'* and extend the formula- 
tion to include forward-flight computations. This 
new formulation yields efficient numerical predic- 
tion of HSI noise without resorting to unnecessary 
or ad hoc simplifications of the FW-H quadrupole 
source term. The mathematical manipulations used 
in this approach are rigorous and depend only on 
the far-field assumption, without approximation of 
the source strength. Numerical time differentiation 
of integrals is avoided throught the development of 
an alternate formulation in which the time differen- 
tiation is done analytically. Preliminary calculations 
with this new formulation demonstrate the potential 
for efficiency and robustness. 

The acoustic analogy approach was chosen be- 
cause of the substantial knowledge base gained in the 
development and utilization of thickness and load- 
ing noise predictions, ba.sed on the FW-H equation. 
Further, the fundamental far-field assumption, 
which is described in the next section, leads to in- 
tegrals of precisely the same form as current thick- 
ness and loading noise calculations; hence, the exist- 
ing numerical algorithms can be used directly. Fi- 
nally, the identification of individual noise compo- 
nents is a unique advantage of the acoustic analogy 
approach. This new formulation has been coded and 
is described in the remainder of this paper. The nu- 
merical results are compared with experimental data 
for both hover and forward-flight conditions. 
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Formulation Development 

The FW-H equation is the most general form of 
the Lighthill acoustic analogy and is appropriate for 
predicting the noise generated by the complex mo- 
tion of helicopter rotors. In differential form, the 
FW-H equation is given by the following inhomoge- 
neous wave equation: 




+ 


{Mf)] - -^[usu)\ 


dxidxj 


[Tamf)] 


( 1 ) 


where p'{x,t) is the acoustic pressure and the three 
source terms on the right-hand side are known as 
the thickness, loading, and quadmpole source terms, 
respectively. The rotor-blade surface is defined by 
the equation / = 0. Note in equation (1) that 
the thickness and loading source terms are surface 
distributions of sources (indicated by the presence 
of the Dirac delta function S(f)). Also note that 
the quadrupole source is a volume distribution of 
sources (indicated by the Heaviside function 
The FW-H equation is valid in the entire unbounded 
space; hence, a formal solution may be obtained by 
using the free-space Green’s function S{g)/ATvr. In 
this paper, we are primarily concerned with the con- 
tribution of the quadrupole source term. 

Farassat and Brentner^^ have shown that after 
some manipulation of the formal solution the noise 
contribution from the quadrupole may be expressed 
as 


47rpg(f,/) = 


c dP 


dU dr 


II 

-oo />0 
-Oo />0 


+ C 


-Tu 


dQ dr 


— oo />0 


where Pg(x, f) is the acoustic pressure due to the 
quadrupole source. The quantity Trr is the double 
contraction Tijvirj, and Vi are the components of 
the unit vector in the radiation direction. In addi- 
tion, dQ is an element of the surface g = 0, which is 
known as the collapsing sphere; hence, equation (2) 
is known as a collapsing-sphere formulation. 


Formulation Q1 

Equation (2) is the starting point for deriving the 
formulation developed by Brentner and Holland. 



Figure 1. Schematic of collapsing sphere that inter- 
sects rotor blade. 


To better illustrate the far-held approximation used 
in reference 14, we hrst give a geometric interpre- 
tation of the collapsing-sphere formulation. The 
collapsing sphere is dehned by the equation g = 
T — t + r/c = 0, where r and t are the source and 
observer times, respectively, and r is the distance 
between the observer position x and the source posi- 
tion y. Because the observer time t and the observer 
position X are held hxed during the integration, the 
solution to 5 = 0 can be interpreted as a sphere cen- 
tered on the observer x of radius r = |x — pj, which 
reduces, or collapses, as r approaches t. A schematic 
is shown in hgure 1. 

An integration over the entire collapsing-sphere 
surface is not necessary in the integrals of equation 
(2) because the Lighthill stress tensor Tij vanishes 
away from the source region. For an observer in the 
far field, the collapsing sphere can be locally approx- 
imated by a cylinder, as shown in figure 2, because 
the source strength reduces quickly in the direction 
normal to the rotor blade. A further simplification 
in acoustic calculations occurs if the observer is as- 
sumed to be in the rotor plane (precisely where HSI 
noise has maximum directivity) because in this case 
integration in the direction normal to the rotor plane 
can be done independently of the observer position. 
Yu et al.'^ were the first to use this far-field approx- 
imation for the evaluation of helicopter rotor HSI 
noise; however, the present work differs in that no 
additional approximation is made to the quadrupole 
source strength. 

The integration over the approximate collapsing- 
sphere surface is carried out in two stages. First, 
integration in the direction normal to the rotor disk 
is accomplished. This step allows us to define the 


3 





Figure 2. Schematic of approximation to collapsing 
sphere in source region. 


quadrupole source strength on the rotor disk as 


Qij — 


/ 


Tij dz 


(3) 


where ^ is understood to be in the direction nor- 
mal to the rotor disk and the z integration is only 
done outside the rotor blade. The new tensor Qij is 
distributed on the rotor disk plane and is expected 
to vanish sufficiently far ahead of the leading edge, 
behind the trailing edge, and off the blade tip. By 
using relation (3), equation (2) can now be written 
as 

/+=0 
g = 0 


+ 


d f f 3Qrr Qii j 

at] J 

/+=0 

g=0 


+ =/ / 
/+=0 
5=0 


where /+ = 0 represents the rotor disk plane, in- 
cluding the blade surface. The intersection of the 
collapsing sphere with the rotor plane results in a 
curve for which we use the notation F. (See figure 
1-) 

The integrals in equation (4) are in the same form 
as the collapsing-sphere formulation for thickness 


and loading noise developed by Farassat.^® For this 

1 

reason, we can apply the relation 


c dr dr _ dS 
sin 0 1 1 — Mr I 


(5) 


to transform equation (4) from a collapsing-sphere 
formulation to a retarded-time formulation. In equa- 
tion (5), dS is an element of the source surface, and 
6 is the angle between the source surface normal and 
the radiation direction. When the observer is in the 
rotor plane, sin^ = 1. This transformation is exact. 

The retarded-time formulation that results from 
the application of equation (5) to equation (4), first 
presented in reference 14, is 


4Trp'Q{x,t) 


+ 


+ 


1 

f \ Qrr 

dt^ J 

1 ir\l-Mr\^ 

/ + : 

=0 

f 

3Qrr Qii 1 

cdt J 

. 1 1 — Mr 1 J 


/+=o 


dS 

1 ret 

dS 

ret 


f r 3Qrr Q»»1 

J ir^\l-Mr\\ 

/ + =0 


dS 

ret 


( 6 ) 


where Mr is the local Mach number of the source 
in the radiation direction and the subscript ret indi- 
cates that the integrand is evaluated at the retarded 
time t — r/c. For convenience, we refer to this equa- 
tion as formulation Ql. (This notation parallels that 
used by Farassat for the thickness and loading for- 
mulation 1.) Although formulation Ql is strictly 
valid only in the rotor plane where sin0 = 1, this 
formulation is expected to provide useful results for 
observers located at small angles above or below the 
rotor plane. For larger angles, the quadrupole con- 
tribution is not expected to be significant. In addi- 
tion, because the far- field approximation only relies 
on the source strength vanishing in the direction nor- 
mal to the rotor, the observer can be relatively close 
to the rotor (e.g., within 1 to 2 rotor radii). 

Evaluation of the quadrupole noise can now be 
completed with substantially less computational ef- 
fort than a direct numerical evaluation of either 
equation (2) or its retarded-time counterpart. Be- 
cause the integrals in equation (6) are of the same 
form as thickness and loading noise, formulation 
Ql can be easily incorporated into existing rotor- 
noise prediction codes such as WOPWOP.^^ Indeed, 
this formulation has been incorporated into WOP- 
WOP by Brentner and Holland. Furthermore, the 
numerical insight, experience, robustness, and effi- 
ciency that has been realized in previous research 
can now be applied directly to the approximate 
quadrupole prediction. 


4 



Formulation QIA 

One drawback of equation (6) is that numerical 
time differentiation of the first two integrals is re- 
quired. If the observer is stationary, then this re- 
quirement is not a problem because the time history 
of the integrals can be easily differentiated numeri- 
cally. If the observer is moving with respect to the 
fluid, as in the case of a wind-tunnel test, the situa- 
tion becomes more complicated because the formu- 
lation requires the observer to be stationary during 
the evaluation of the integrals. Predictions with a 
moving observer are possible by adjusting the ob- 
server position at each time in the acoustic-pressure 
time history; however, three evaluations of the inte- 
grals are needed to perform a second-order central- 
difference approximation to the time derivatives at 
each observer time. These extra integral evaluations 
are unnecessary if the time derivatives are taken in- 
side the integrals analytically. 

Although the derivation of a formulation with the 
time derivatives inside the integrals is not difficult, 
it is quite tedious. As a starting point, the integrals 
in equation (6) can be labeled conveniently as 7 i, 
l2, and /a, respectively. Thus, equation (6) can be 
written as 

47 Tp'Q{x,t) = ~(li) + ^{h) + h ( 7 ) 

Now, to take the observer time derivatives inside the 
integrals we use the relationship given by Farassat 

1 Q 

and Sued: 


d 


1 d 


dt 

X 

1 — Mr dr 

X 


This result is obtained by using the definition g = 
T — t + vjc = 0 , given the fact that r is a function of 
source time r through the source position variable 
y. Two applications of relation (8) on I\ yield 


/ + =0 



/+=0 


where A'l, /l2, and A'a are defined as 

o/ dM. \2 

r.' _ dT^ , dT > 

^ - MrY c^r{\ - Mrf 

cr2(l - MrY r3(l - Mrf 


3 ^ , 

c'^r{l-MrY cr\\-Mrf 
1 

c2r(l - Mrf 
The time differentiation taken inside I-2, yields 

= J ( 10 ) 

/ +=0 



I<2 

and 

I <3 


with A4 and A's defined as 

^ 2Mr 

cr^{l-MrY r^{l-MrY 

1 

cr2(l — MrY 

Terms K\ through A's depend only on the kinemat- 
ics of the source motion and have been determined 
by using relation (8) together with 


AT 4 = 

and 

As = 


5 ; = -CM, ( 11 ) 

The values of A'l through A's can be evaluated by 
utilizing the definitions of the source-time deriva- 
tives of Mr, which are written as 


dMr 

dr 

= Mr + -{M^ + M^) 
r 

(12) 

d'^Mr 

3 c • ■ 

( 13 ) 

dr^ 

= Mr + —{Mr - MiMi) 
r 


3 c Mr , , ,2 »r2\2 

+ — ~{M^ - M^Y 



These relations follow directly both from the fact 
that dy/dr = v = cM and from the definition of r. 

The derivation would be complete at this point 
if the time-dependent values of Qrr and Qa were 
given as input for each point in the source region. 
To compute Qrr, however, requires knowledge of the 
observer position, which is unrelated to the flow- 
field calculations. For this reason and to take full 
advantage of the in-plane observer assumption, we 
assume that the input data do not require knowledge 
of the observer position. Hence, the determination of 
the source-time derivatives of Qrr and Qu remains. 
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By recalling that 


Qrr — Qij'^i^j 


and 


drj _ c{fjMr - Mj) 

dr r 

we can easily to show that 

dQrr A I ‘2>c(^MrQrr Qmm') 

5 — Vrr H 

OT r 


(14) 


d^Q„ _ 

dr^ 


(15) 


i&C • ‘ * 

+ ~{MrQrr + 2MrQrr ~ SQjVfr ~ Q Mr) 

2c^ 

+ -r[QMM + - M^)Qrr ~ 4M,QMr] 


and, by definition, 

^ = Qii (16) 

The notation used here is defined in the Notation 
section; however, the reader should understand that 
the dot (■) implies source time differentiation on the 
variables shown with the dot. A dot on the main 
variable does not imply differentiation of any of the 
associated vectors implied by the subscripts. Sub- 
scripts other than i and j are a shorthand for the 
inner product of the main quantity with the vector 
represented by the subscript. 

The final form, which is obtained by combining 
the previous results, becomes 


i^p'Q(x,t)= I + dS (17) 

^ J Lc-^r cr-^ H Jre< 


/ + =0 


where 

Qrr ^^rQrr T SAff-^rr 

(1 - Mrf ^ (1 - MrY 

2M^Qrr 
(l-Mr)^ 

j. _ —Qii + "^QMr ^rQii 

“ (1 - MrY ~ (l-Mr)^ 

3[(1 - M^)Qrr - 2MrQMr “ MiMiQrr] 

6 Mr(l - M^)Qrr 

(1 - 


„ 2Qmm - (1 - M^)Qii 6(1 - M^)Qur 

3(1 - M^Q„ 

(1 - Mrf 

Equation (17), together with the definitions of i^ri> 
Kr 2 , and Krs, will be referred to as formulation QIA 
and is the main result of this paper. The designation 
QIA is intended to parallel that of Farassat’s formu- 

17 .. 

lation lA, in which the observer time derivative is 
taken analytically inside the thickness and loading 
integrals. Formulation QIA does not require nu- 
merical time differentiation of the integrals and, as 
a retarded-time formulation, is well suited for sub- 
sonic source motion. Aside from the problem geome- 
try, only the time-dependent value of Qij is required 
as input. 

Note that the definition of the Lighthill stress ten- 
sor Tij, as well as that of Qij , requires that the fluid 
velocities be specified in an inertial reference frame 
fixed to the undisturbed medium. In this frame, the 
fluid velocities are perturbations about the station- 
ary medium. Often, CFD solutions are computed in 
a body-fixed reference frame; hence, the motion of 
the body must be removed. Finally, remember that 
the vector M (and all related quantities such as Mr) 
is the source velocity rather than the fluid velocity. 

The remainder of this paper will discuss the nu- 
merical implementation of equation (17) and present 
some initial predictions. 

Numerical Implementation 

The numerical calculation of the quadrupole noise 
has been divided into two stages: a preprocessing 
stage in which the integration of the Lighthill stress 
tensor in the normal direction, indicated in equa- 
tion (3), is carried out, and an evaluation stage in 
which the quadrupole contribution to the acoustic 
pressure specified in equation (17) is determined. 
Both the preprocessor and the acoustic calculation 
are described in this section. More information on 
the program details can be found in reference 14. 

Preprocessor 

Although the evaluation of Qij can be performed 
independently of the observer position and the re- 
tarded time, the preprocessor must read in the CFD 
solution, interpolate the solution at the necessary 
quadrature locations, and then perform the numer- 
ical quadrature in the direction normal to the rotor 
disk. The preprocessor needs knowledge of both the 
CFD grid topology and the solution format. In the 
initial implementation used for this work, the inter- 
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polation of the CFD data is two dimensional and 
is done one radial station at a time. For a given 
radial station, data are interpolated to quadrature 
points needed for composite Gauss-Legendre inte- 
gration, on lines normal to the rotor plane. The 
lines are uniformly distributed in the chordwise di- 
rection. A two-dimensional linear least-squares in- 
terpolation is used to interpolate the density, mo- 
mentum, and energy at each quadrature point. The 
Lighthill stress tensor T{j is evaluated with the in- 
terpolated data. The value of Qij on /"'■ = 0 is 
determined at each chordwise location before mov- 
ing to the next radial station. The entire process is 
repeated for each source time and is stored as input 
to the acoustic calculation stage. 

Noise Calculation 

The quadrupole noise calculation indicated in 
equation (17) is carried out in a modified version of 
WOPWOP,^^ called WOPWOP-I-. The quadrupole 
noise calculation is distinct from but parallels that of 
thickness and loading noise in the WOPWOP code, 
so the discussion here will be brief. The WOP- 
WOP-f code used in this paper is essentially the 
same as that used in reference 14, with two ex- 
ceptions. First, the code was modified to use un- 
steady quadrupole source loading. Second, formula- 
tion QIA was implemented to enable forward-flight 
wind-tunnel predictions. 

One unique aspect of the quadrupole noise calcu- 
lation is that the integration surface extends ahead 
of the leading edge, behind the trailing edge, and 
off the rotor-blade tip. Because of this difference in 
the integration surface, the quadrupole noise calcu- 
lation is a separate subroutine in WOPWOP-I-. The 
panels that extend off the tip can be a problem for 
a retarded-time formulation because rotating pan- 
els move supersonically if the radius of rotation is 
great enough. Rather than implement specialized 
code logic to handle the singular integrals caused by 
supersonic panel motion, only panels with subsonic 
speeds are included in the present calculations. 

Because the CFD solutions used in this work are 
output in a coordinate system that is always aligned 
with the rotor blade, coordinate transformations are 
used to transform Qij into the ground fixed frame 
used for all acoustic calculations. The numerical in- 
tegration over each panel is carried out for all ob- 
server times before proceeding to the next source 
panel to reduce the retarded-time differences be- 
tween subsequent evaluations. Unlike the thickness 
and loading algorithm in WOPWOP, which allows 
the user to select the number of computation panels, 
the quadrupole calculation in WOPWOP-f uses the 


computational grid from the preprocessor without 
modification. Gauss-Legendre quadrature is used in 
both the chordwise and radial directions for the nu- 
merical integration over each panel. The retarded 
time is determined by using the modified false posi- 
tion method. 

An adaptive numerical integration scheme is used 
to compute the contribution of the quadrupole noise 
from each panel. The motivation for using an adap- 
tive scheme is that, although the quadrupole source 
strength may vary essentially linearly over a panel, 
experience has shown that when the retarded time 
is properly taken into account a larger number of 
quadrature points is needed if the effective area (i.e., 
the emission area AE, also called the acoustic plan- 
form) of the panel is large. The panel area, panel ve- 
locity, and observer position all determine the acous- 
tic planform of the panel. When the acoustic plan- 
form of the panel is large, more quadrature points 
are expected to be needed to perform the integra- 
tion over the panel because the integrated function 
may have more variation. An adaptation parameter 
given by 

where Mr is the Mach number of the panel in the 
radiation direction and A is the physical area of the 
panel, is used to determine the number of quadrature 
points for a particular panel; N is scaled such that 
1 < TV < 20. As Mr approaches unity or if the 
area of the panel is large, the adaptation parameter 
becomes large. Similarly, if Mr is small or the area 
of the panel is small, then the adaptation parameter 
becomes small. For this paper, only the number of 
points in the chordwise direction has been varied 
because little difference is found when more points 
in the radial direction are used. 

Sample Predictions 

In this section, two representative calculations are 
performed to demonstrate the new formulation and 
to provide some indication of the efficiency and ro- 
bustness of the new quadrupole prediction method. 
Predicted acoustic-pressure time histories will be 
compared with meeisured data for both hover and 
forward-flight conditions. 

Hover 

A model scale rotor test, conducted by Boxwell 
et al.^^ in 1978 and later repeated by Purcell^^ 
in 1988, was selected for validation of the present 
theory and code. Both experiments were designed 
specifically to measure HSI noise generated by a non- 
lifting helicopter rotor in hover. The rotor was a 
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Figure 3. UH-IH model rotor operating in hover 
with Mh = 0.88. ■ experimental data; pre- 
dicted acoustic pressure; quadrupole compo- 

nent of predicted acoustic pressure. 

one-seventh scale model of a UH-IH main rotor with 
straight, untwisted blades. The model rotor utilized 
an NACA 0012 airfoil section. The rotor radius R 
was 1.045 m with a chord of 7.62 cm. The model 
was run at several high-speed hover conditions with 
low thrust. Comparisons of measured and predicted 
acoustic pressure are made for an in-plane micro- 
phone located 3.09/2 from the rotor tip. 

For the hover noise calculation, an Euler solution 
reported by Baeder et al.“^ was used as input. The 
Euler calculations were performed on a C-H grid; 
only the lower half of the grid was used in the calcu- 
lations by taking advantage of the symmetry of the 
problem. The Euler calculations required approxi- 
mately 80 min of CPU time on a Cray Y-MP. Details 
of the Euler calculations can be found in references 
21 and 22. 

In figure 3, the predicted quadrupole component 
and the total acoustic pressure are compared to the 
experimental data for a tip Mach number Mh = 
0.88. For this case, the flow around the UH-IH ro- 
tor blade is transonic but not delocalized; neverthe- 


less, the quadrupole contribution to the total sig- 
nal is substantial. The good agreement between the 
current prediction and the experimental data is rep- 
resentative of HSI noise predictions with this new 
quadrupole formulation when the rotor speed is be- 
low the delocalization Mach number. Although the 
result here is a new computation using formulation 
QIA, this result is indistinguishable from the pre- 
vious calculation by Brentner and Holland, which 
used formulation Ql. For this preliminary calcula- 
tion, approximately 6 sec of CPU time on a scientific 
workstation was required for the preprocessor execu- 
tion, and an additional 40 sec was needed to run the 
combined thickness, loading, and quadrupole noise 
prediction. 

Forward Flight 

To demonstrate the forward-flight capability of 
the new formulation and code, a comparison of pre- 
dicted and measured results for a four-blade swept- 
tip rotor tested in the Duits-Nederslandse Windtun- 
nel (DNW) was made. For this comparison, a micro- 
phone located in the rotor plane at a rotor azimuth 
of Ip = 150° was utilized. The model rotor was oper- 
ated at an advance ratio /z = 0.32 and an advancing 
tip Mach number Mat = 0.933. The experiment is 
described in the report by Visintainer et al.^^ 

For this prediction, the full potential solver 
FPRBVI^^ was used to compute the unsteady flow 
field around the rotor. For this developmental cal- 
culation, the FPRBVI code was modified to store a 
solution file at every degree of rotor azimuth. These 
files were then read by the preprocessor program to 
compute the time-dependent values of Qij . Manipu- 
lation of the three-dimensional time-dependent data 
was computationally demanding in terms of CPU 
time and storage. For this case, approximately 1100 
CPU sec on a scientific workstation was required by 
the preprocessor. Both computer time (by reducing 
input and output requirements) and storage could 
be reduced significantly if the preprocessor was in- 
cluded as a data post-processing subroutine in the 
CFD code. 

The results of the forward-flight prediction are 
shown in figure 4. Again, the experimental data are 
compared with the predicted acoustic pressure; the 
quadrupole contribution is also shown to indicate its 
relative magnitude. Although the CFD calculation 
used a rather coarse grid, the agreement is good. 
Further work is necessary, however, to determine 
whether the CFD solution is fully converged and to 
determine the sensitivity of the quadrupole predic- 
tion to various input parameters. This develop men- 




Figure 4. Contemporary design, four-blade model 
rotor operating in forward flight; ^ — 0.32 and 
Mat = 0.933. ■ experimental data; pre- 
dicted acoustic pressure; quadrupole compo- 

nent of predicted acoustic pressure. 

tal calculation required approximately 190 CPU sec 
on a scientific workstation. 

Concluding Remarks 

The new far-field quadrupole noise formulation 
presented in this paper provides an efficient method 
for computing high-speed impulsive (HSI) noise. 
The method is efficient in the sense that the vol- 
ume integration normally associated with the Ffowcs 
Williams-Hawkings (FW-H) quadrupole has been 
reduced to a surface integration of the same form as 
thickness and loading noise. Taking the time deriva- 
tives inside the integrals also improves the efficiency 
of the numerical method by avoiding extra integral 
evaluations required by numerical time differentia- 
tion. Even though the code has not been optimized, 
the acoustic predictions are already relatively fast 
for both hover and forward-flight predictions. 

The new formulation should also prove to be ro- 
bust. Brentner and Holland^^ have demonstrated 
that the far-field approximation gives the correct sig- 
nal amplitude, even for high-speed cases in which the 


present retarded-time formulation does not strictly 
apply. In these cases, supersonic quadrupole sources 
are ignored in WOPWOP-f ; the signal shape suffers 
and is incorrect. This graceful failure is character- 
istic of the robust nature of retarded-time formu- 
lations. The efficiency and robustness of this new 
quadrupole formulation make the acoustic analogy 
an attractive option for the prediction of HSI noise 
and will help to maintain the usefulness of rotor 
noise prediction codes based on the FW-H equation. 
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